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ABSTRACT 


In this dissertation we have studied the balaruced netxj>ork 
flow problem CBNFPD . In the balanced network flow problem, the 
objective is to balance the weighted arc flows in a distribution 
network, by minimizing the difference between the largest and the 
smallest weighted arc flow. This type of problem arises in 
military, industrial, and economic applications requiring equitable 
distribution of resources. 

The balanced network flow problem can be transformed to 
the standard form of the linear programming problem by introducing 
2m+2 additional variables and 2m additional constraints, where m is 
the number of arcs in the network. This transformation is not 
desirable from computational point of view for large size of m as 
it increases the size of problem substantially. In this thesis, we 
present two algorithms, namely. parametric algorithm, and binary 
search algorithm for the balanced network flow problem, which do 
not require enlargement of the problem. The parametric algorithm 
utilizes the relationship between BNFP and mxnimxxx network flow 
problem CMNFP^ and parametrically solves the latter problem. The 
binary search algorithm uses binary search technique to reach near 
the optimal solution. Both the algorithms are coded in FORTRAN-77 
and tested extensively with random networks. We compare the 
computational performance of these algorithms with that of the 
simplex algorithm, applied to the enlarged linear programming 
problem. Our computational investigations show that the parametric 
algorithm is computationally superior than the other two 


methods . 



CHAPTER I 


INTRODUCTION 


1.1 INTRODUCTION 

I*t IS very impor'tan't -to distributee a given amount of 
resources &<^VLitably in many situations. By table distributzorif 

we mean a distribution scheme, which allocates a resource to 
different centers proportional to their demands. Obviously, 
equitable distribution arises, whenever there is a scarce resource 
or nature of the distribution network is such that, it is not 
possible to satisfy exact requirements of demand centers. Under 
such conditions, one may like to distribute the resource in such a 
fashion, that each center get a share proportional to their 
requirement Cequi table distribution!). However, the distribution 
network may not allow the perfect equitableness, when capacity 
constraints are present. So a distribution scheme might attempt 
to, 

CID maximize the smallest proportion; or 

C2D minimize the largest proportion; or 

C3D maximize the smallest proportion and minimize the 
largest proportion simultaineously. 

In this thesis, we discuss a distribution scheme which 
tries to attain the equitableness by minimizing the difference 
between the largest and smallest proportions. 



1.2 


PROBLEM DESCRIPTION 


Lot, G = CN,A3 bo a dirocted nobwork wibh cost c. lower 

i ^ 

bound 1. . and capacit-y u. . associated with each arc Ci,jD e A. We 
associate with each node i € N> an integer number bCiD , represent- 
ing its supply or demand. If faCiD > O, then node i is a sxxpply 
Tuodba^ if bCiD < O, then it is a demand node, and if bCiD =0, then 
node i is a transshipm&nt nod^. Let n = |N| and m = |A[. Let A ^ 

A. Then we define the B<xlainc&d N&iwoi^h Flow Prablom CBMFF^ as 
follows : 


Minimize -( 


Max. i c. .X. . I ~ Mi 
* 

Ci , jO € A Ci , j:> 


.n. -I c. .X. . I 

I J 


sub j ect to 

J x.j - 2 = bCi3, V i ^ N 

< j : Ci,jD € A > < j : Cj.iD € A > 

1 . . < X. , < u. , , V Ci,jD € A. 
ij ij iJ 


Ci. iy 


ci.2:> 


Ci.SD 


We refer to the vector x = Cx^^D as the flow in the 
network. Each constraint in Ci.2^ implies that the total flow out 
of a node minus total flow -into that node must equal to the net 
supply / demand of the node. We henceforth refer to these constra- 
ints as the mass balance const Taints. The flow must also satisfy 
the lower bound and upper bound constraints Ci.3^, which we refer 
to as the capacity constraints. 



1.3 


ASSUMPTIONS 


Throughout- t-he material in this thesis, the following 
assumptions have been made. 


CID 


C2D 


Although, in general, arc cost ‘”j_j can be in the 
interval C— oo, coD we have considered the common 
practical case of 2 0 for each Ci,JD e A. 

To achieve simplifications in the presentation we 


assume A = A . However we note that this assumption 
is only for ease of presentation and later we 
explains how this assumption can be relaxed. 


We assume in the rest of the thesis that these assumptio- 
ns are satisfied unless we state it otherwise. 


1 . 4 - APPLICATIONS 

Resource sharing is a problem arising in numerous 
application settings. Many important problems such as distributing 
relief to famine stricken region* allocation of federal govern- 
ment’s funds to provinces, distribution of essential commodities in 
retails, requires sharing of resources. This section illustrates 
two scenarios where balanced network flow problems may arise in 
practice. 


COAL STRIKE PROBLEM [13 


Consider the situation during a prolonged coal strike. 
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Some nonunion mines could be producing during the strike. Let Q be 
the set of nonunion mines and a^ be the amount of coal produced by 
the mine q for each q e Q. The problem is to distribute the 
limited coal supply produced by the nonunion mines equitably among 
the power companies. Distributing each power company the same 
amount will ignore the size of the company and would allocate 
smaller companies more coal while the larger one will face short- 
age. Thus, the distribution scheme must consider the relative size 
of the power companies. Let P denote the set of the power compan- 
ies and d^ for p e P represent the amount of coal used by the power 
company p during normal times. Let fCpD for p e P denote the 
amount of coal the distribution scheme assigns to power company p. 
So fCpDx^d^ is the proportion of the amount of normal time demand 
supplied to power company p. An equitable distribution scheme 
would attempt to equalize the proportion. Perfect equalization may 
not be possible due to capacity constraints of the distribution 
network. So one can try to balance the proportion by minimizing 
the difference between the largest and smallest proportion. The 
problem can be modeled as a BHFP by appropriately defining the 
network and the terms. Consider a distribution network connecting 
the nonunion mines and power company as in Figure 1.1. Let s be a 
super source node and t be a sup>er sink node. Join s to node q by 
an arc Cs,qD for each q e Q. Similarly, join p to node t by. an arc 

Cp.tD for each p e P. Now t-he disbribut-ion scheme can be sbabed in 

# 

the form of equabion Ci.Syt Ci . 3J> where A = -C Cp, bD : p 

P > , c , =1 /d , suppl y ab node s i s bC sD = Y a , demand ab b i s 

pb p Z# q 

bCbD = -bCsD and bCi3 = O for obher nodes. 
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Notwork Mocoi of Coal Striko Example 


Figure 1 .1 



FLOOD RELIEF PROBLEM 


Consider a region affected by the flood. There are some 
relief shelters in that region where the people of the affected 
areas have taken shelter. ' Let us denote the shelters as r^ , r^, 

. . . ,r . Suppose there are m relief agencies a. , a^, . . . , a 

operating in these areas. Since each relief agency has limited 
supply, for effectiveness of their mission they would like each 
person of all the shelter to get an equal share of the supply. It 
may not be possible to supply each shelter an amount of supply 
proportion to the number of people staying in it because of the 
nature of transportation network. This problem is similar to the 
coal strike problem if we think of the shelters as the power 
companies and relief agencies as the nonunion mines. 

1.5 MOTIVATION OF THE THESIS 

The BNFP can be formulated as a Linear Programming 
Problem CLPD as below ; 


Minimize Cz - yD 


Ci . 4:> 

subject to 



y X. . - 

Z XJ 

Tx.. =bCiD, Vi eN, 

A 

Ci.sy 

< J :Ci,jD € A > 

<j:Cj,i5€A> 


C.jXij <2, V Ci, 

j:> € A, 

ci.ey 


j:> € A. 

ci.ry 

l.j < X.J < U.J. 

V Ci, jD € A. 

Ci.3:> 


This linear programming formulation can be converted to 



■the st.andard form by in'Lroducing 2m slack variables. The enlarged 
LP can be solved by any of bhe existing simplex based algorithms. 
These algorithms however consider the constraints Ci , 6:> and Ci.7> 
explicitly which results xn the following disadvantages. 

Cl 3 As the value of m increases, the size of the basis 
increases resulting in large storage requirement and more computat- 
ion per iteration. 

C2D The basic solutions are highly degenerate, specially 
in the first few iterations. So the algorithm may go on iterating 
without increasing the objective function. 

C3]) The coefficient matrix of the mass balance constrain- 
ts possesses special structure. Only 2m out of its nm entries are 
nonzero. All of its nonzero entries are +i or -1 and each column 
has exactly one +1 and one -1. However, the simplex based algorit- 
hms for LP does not exploit its special structure to achieve 
theoretical and computational simplifications. 

In this dissertation, we present a simplex based algo 
rithm for the BNFP which takes care of all the disadvantages stated 
before. It considers the additional constraints Cl , 6S> and Cf. 7J> 
implicitly and exploits the special structure of the constraint 
matrix. The algorithm is based on the relationship of BNFP with 
the Minimax Network Flow Problem CMNFPD and parametrically solving 
the latter problem. 

16 NOTATION AND DEFINITIONS 


In this section w© state without proof some of the basic 



properties of graphs. The notation and conventions which are used 
in the thesis are also given for the sake of completeness. 


A directed graph G = CN,AD consists of a set, N, of 

nodes, and a set. A, of arcs whose elements are ordered pairs of 

distinct nodes. A directed network is a directed graph with 

numerical values attached to its nodes and/or arcs. As before, we 

let n = |N| and m = |M|. We associate with each arc Ci,jD e A a 

cost c. a capacity u. and a lower bound 1, .. We assume through— 
ij 

out that l^j > O for each Ci,JD e A. Frequently we distinguish two 
special nodes in a graph : the source s and the sink t. 

An arc Ci,j3 has two end points, i and j. The arc Ci.jD 

is incident to nodes i and j. We refer to node i as the tail node 

and node j as the head of arc Ci,jD, and say that the arc Ci,j:) 

emanates from node i. The arc Ci,jD is an outgoing arc of node i 

and an incoming arc of node j. The arc adjacency list of node i, 
ACiD , is defined as the set of arcs emanating from node i, i . e. , 
ACiD = <Ci,jD e A : j e N >. The degree of a node is the number of 
incoming and outgoing arc incident to that node. 

A directed path G = CN,A5 is a sequence of distinct nodes 


and arcs i^, i^, ' '^^r-1 ’^r^ ’ ^r satisfying the 

property that Ci , i e A for each k = 1,. . , i 1. An undirected 

path is defined similarly except that for any two consecutive nodes 
i^ and the path, the path contains either arc 

arc Ci, ,i, 3. We refer to the nodes i^, i^,. . »i__^ a-s the internal 

nodes of the path. A directed cycle is a directed path together 

with the arc Ci ,i^3 and an undirected cycle is an undirected path 

r 1 

together with the arc Ci^,i^I> or Ci^^.i^D. 

We often use the terminology path to designate either a 



directed or an undirected path, whichever is appropriate from the 
context. If any ambiguity arises, we explicitly state directed or 
undirected path. For the simplicity of notation, we often refer to 
a path as the sequence of nodes i^-ig—. . -i^ when its arcs are 
apparent from the problem context. Alternatively, we sometimes 
refer to a path as a set of C sequence of !> arcs without mention of 
the nodes. We use similar conventions for representing cycles. 

# # # i 

A graph G = CN,AZ)isa subgraph of G = CN,AZ) if N^A 

* * 9 i 

and A ^ A. A graph G = CN,ADisa spanning subgraph of G = 
CN,AD if n' = N and a' ^ A. 

Two nodes i and j are said to be connected if the graph 
contains at least one undirected path from i to j. A graph is said 
to be connected if all pairs of its nodes are connected : other- 
wise, it is disconnected. We always assume graph G is connected 

unless and otherwise stated and hence m > n-1 . We refer to any set 

/ 

Q S A with the property that the graph G = CM,A-QD is disconnect- 
ed, and no subset Q has this property, as a cutset of G. A cutset 
partition the graph into two sets of nodes, X and N— X. We shall 
alternatively represent the cutset Q as the node partition CX,N— X3. 

A graph is acyclic if it contains no cycle. A tree is a 
connected acyclic graph. A subtree of a tree T is a connected 
subgraph of G. Arcs belonging to T are called nontree arcs. A 
spanning tree of G = CN,AD has exactly n— 1 tree arcs. A node in a 
tree with degree equal to one is called a leaf node. Each tree tias 
at least two leaf nodes. 

A spanning tree contains a unique path between any two 
nodes. The addition of any nontree arc to a spanning tree creates 
exactly one cycle. Removing any arc in this cycle again creates a 
spanning tree. Removing any tree arc from a spanning tree creates 



"two sub'troes. Arcs whos© ©nd poiii'ts b©long *to bh© bwo dif'f©r©nt- 
subtrees constitute a cutset. If any arc belonging to this cutset 
is added to the subtrees, the resulting graph is again a spanning 
tree. 


1.7 ORGANIZATION OF THE THESIS 


bel ow. 


A brief chapterwis© description of the thesis is given 


In Chapter 2, the existing literature related to the 

problem is discussed. The literature is mostly in the area of CID 
Maxi min linear programming problem, C23 Mini max linear programming 
problem, C3D Mini max transportation problem, C4I> Mini max network 
flow problem C53 Bottleneck transportation problem, C63 Bottleneck 
linear programming problem, C73 Balanced linear programming 
probl em. 

In Chapter 3, we describe in detail the parametric 

algorithm and binary search algorithm to solve the balanced network 
flow problem. The basic ideas and concepts of these algorithms are 
described, followed by their stepwise descriptions. We illustrate 
the steps of these algorithms with a numerical example. 

In Chapter 4, we give a detail description of our 

implementations of these algorithms. We also report the computati- 
onal investigations of these algorithms and the simplex algorithm. 
The results of our computational testing for different sizes of 

random networks are also presented. 

We give the list of references and bibliography at the 


end. 



CHAPTER II 


LITERATURE SURVEY 

In this chapter, a brief literature survey in the area of 
minifttax and maximin optimization is presented. Over the years, 
many researchers have studied the minimax and maximin optimization 
problems and their results have been reported in the literature. 
We will go through some of these results in the following discuss- 
ion . 

2.1 MAXIMIN LINEAR PROGRAMMING PROBLEM 

Maximin Linear Programming Problem has been studied by 
Kaplan C2], Brown C33, Gupta and Arora C4], Wu and Posner C3]. A 
Maximin Linear Programming Problem is described as an optimization 
problem with a maximin objective function subject to the linear 
constraints. Mathematically, it can be represented as follows : 

Maximize Z = Minimum Cc x } 

1 1 

l:l€(l....n> 

Subject to 

A X = b. 

X > 0. 

where c^ and x^ are scalars for i = l,..,n, X is a column vector of 
Xj, A is an mxn matrix and b is a column vector of m rows with 
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elements b . 

1 

Although Maximin Linear Program is mathematically not 
difficult and is similar to other nonlinear problems that can be 
solved by linear programming , Kaplan C23 has described a simplifi- 
ed solution procedure when certain constraint conditions exists. 
He has obtained a simple solution for > 0 when the following 
condition holds ; There exists a feasible solution X, such that 
c^Xj = t V 1 = 1, .. ,n. Sucn a solution is called a ray solution 
and if it exists. then it clearly solves the Maximin Linear 
Programming Problem. He has also pointed out many military ana 
economic applications of the problem. 

Brown C53 termed the linear programming problem with 
maximin objective function as *'lin&aT sharing problem.". He has 
Developed an algorithm which alternatively calculates a new global 
upper bound on the objective function value and then determines if 
a feasible solution exists that meets the new objective. 

Gupta and Arora C43 have extended the work by Kaplan C23 
to consider more general cases. Wu and Posner C53 have studied 
several aspects of maximin linear programming problem. -Their 
algorithm can be divided into two parts. The first part finds an 
initial feasible solution to maximin linear programming problem. 
It is essentially an extension of Kaplan’s C23 ray solution and 
whenever such an procedure is not applicable, phase I of the 
simplex method has been used to find an initial feasible solution 
of the set Cx|Ax<b;x>0>. The second part is an iterative 
method for searching the optimal solution. 

Jocobsen C63. Porteus and Yonmark C73 and Brown [83 have 
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ained mi n imi zat i ons . Ahuja Clll has developed two algorithms, one 
parametric algorithm which solves linear programming problem witn 
parametric upper bounds and the primal dual algorithm which solves 
a sequence of related dual feasible linear programming problems. 
Similar to the minimax linear programming problem, the maximin 
linear programming problem can be transformed to a linear program 
by introducing additional constraints. These additional constr- 
aints can be viewed as parametric lower bound constraints and can 
be considered implicitly. In view of this, parametric and primal- 
dual algorithms of Ahuja Cll] can be easily adapted for the maximin 
linear programming problem. 

2.3 MAXIMIN NETWORK FLOW PROBLEM 

Minimax / maximin optimization problem in transportation 
and network context has been studied by various researchers. Brown 
C12] has considered a flow circulation problem with maximin 
objective functions. He has partitioned arcs in the network into 
regular arcs and trade off arcs, where each trade off arc has a 
nondecreasing trade-off functions associated with it. All arcs 
have lower and upper bounds on their flows while the value of the 
smallest trade-off function is maximized. Each cut in the network 
forms a knapsack sharing problem which leads to an optimality 
condition similar to maxfiow mincut theorem. Brown's flow circul- 


ation problem can be stated as ; 
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Maximize i Minimum 4 > [x. .3 1- 

I ^ij ij J 


Subject to 


< i , J ) € T 


I ’^.1 ° I ’‘k.’ ‘ ^ ^ 

j € V k e V 

< x^j < Ujj. V (i.j) eRU T 


where 0 < < u^^ for each arc (i,j) e R u T and V is the set of 
vertices, R is the set of regular arcs and T is the set of trade 
off arcs, X is the flow on arc (i,j>, and u and 1 represent 
upper and lower bounds on arc flows. 


2 . 4 . 


MINIMAX NETWORK FLOW PROBLEM 


A minimax network flow problem on a network G = (N,A> can 
be represented mathematically as follows ; 


z = Minimize s maxi mum <c 
'( 1 , ) > € A 




X >) 
I ] 1 ) J 


subject to 


■ 

( 1 , J ) € A 

y X 

L s 1 
( S , 1 ) € A 


I X,. ^ 

< J , 1 ) € A 




* 

= V , 


( i , t ) e A 


Oix <u , (i,j)€A, 

1 J I J 


where s and t respectively denote the source and sink nodes of the 
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network, ^ ^ , c.^ are flow, capacity and cost (weight) 

associated with the arc <i,j) respectively. 

Stanat and Mago C153 have considered minimax flow in 
linear graphs where a linear graph is defined as a connectea 
acyclic graph with degree of each node either one or two. They 
have considered a flow problem in which a convnodity flows from 
source nodes to sink nodes. Each source node has a specified value 
denoting amount of commodity to be disposed of and each sink node 
has a specified value denoting amount of commodity it can absorb. 
A solution is optimal, if largest flow along any edge is as small 
as possible. They have presented an OCnD algorithm for the problem. 

Ichimori et al . Cl-4] have developed a polynomial time 
algorithm for the minimax network flow problem (MNFP). They have 
called the problem as M>&igh.t&ci minimax flow problem. The polynomi- 
al bound of the algorithm is based on the integral assumption of 
capacity, cost (we i ght ) , and flow value. They have used capacity 

modification technique which treats capacity as a parameter. Using 

* 

the concept of maxflow mincut theorem, an interval on z is 
obtained. They have used binary search technique in that interval 
to obtain z . The time complexity of their algorithm is shown to 

o 

be CXpn D , where n is the number of nodes in the network. p = 

[’log(K>] ,K = Max. (c u ), and [’x'] means the smallest integer y 
(1,J)€A‘^ 

such that y > X. 

Ichimori et al. C153 have extended their work in C143 to 
find polynomial algorithm for MNFP with real valued capacity, cost 
and flow. They have used the capacity modification and strategy 
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for ratio minimization technique of Meggido’s C163. Complexity of 
this algorithm is shown to be OCCCn,iTO^D. where CCn,iTO is the the 
time complexity of maxfiow on a graph with n nodes ana m arcs. 
Currently CCn, nO = OCnm logCn^XirOD Cl 73, 

Finke C183 and Tamir C193 have investigated MNFP with 
integer valued flows. Tamir has demonstrated that bottleneck 
integral flow can be obtained by solving at most k problems with 
real valued variables, where k is the number of discrete variables 
in the problem. Finke C193 has developed a primal algorithm for 
the integer minimax network flow problem. 

2.5 MAXIMUM BALANCED FLOW PROBLEM 


Maximum balanced flow problem has been investigated by 
Minoux C203. Maximum balanced flow problem is a problem of finding 
a maximum flow in a two terminal network such that eacn arc flow 
value is bounded by a fixed proportion of total flow value from 
source to sink. It can be stated mathematically as follows : 

Maximize v 
subject to 

j : <i.j> 

y ^ 

4, SI 

< s , 1 > € A 

0 < X < 

1 ) 


^ x^^ =0, 1 € N - Cs,t}, 

€ A j : < j , i > € A 

2 "u = 

< 1 . t ) « A 


a V , 
1 ] 
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where s and t are source and sink node respectively, a is called 

1 J 

the balancing rate of the arc (i.j). 

Maximum balanced flow problem has been motivated by 
reliability consideration of communication network. If a balanced 
flow IS assigned in a network it is guaranteed that at most the 
fixed proportion of flow value is choked when some arc is broken by 
failure. Fujishige et al . C21] have pointed out the equivalence of 
the maximum balanced flow problem and MNFP. He has shown that MNFP 
can be reduced to a variant of the maximum balanced flow problem. 

2.6 MINIMAX TRANSPORTATION PROBLEM 

t 

Minimax transportation problem <MTP> is a variant of 
the classical transportation problem and can be stateo as follows : 


Minimize Z<x) 

subject to 


Maximum £c x } 
1 J 1 ) 

< 1 , j ) « A 


2 

1 

2 


i j 
€ J 


t 1 


i € I 


2 s 


= a^ . V I € I , 


= b J . V J € J, 



J € J 


where I = C l,...m > is a set of origin nodes. J = C m+1 , . . ,m+n >is 
a set of destination nodes and c^^ is the cost of unit shipment on 
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arc ( 1 , j > . 

Minimax transportation problems have been studied by 
Goldman [Z2.1 and Ahuja C25]!. Goldman C2Z3 has considerea minimax 
transportation problem when = 1 V (i,j> sA . Ahuja C233 has 

developed a parametric algorithm and a primal dual algorithm for 
MTP. The parametric algorithm solves a transportation problem with 
parametric upper bounds and primal dual algorithm solves a sequence 
of maximum flow problem. Worst case complexity of primal-dual 

4 

algorithm is shown to be OCn 3. He has also described a polynomia- 
ily bounded algorithm to find an integer optimum solution. 


2.7 BOTTLENECK TRANSPORTATION PROBLEM 


The bottleneck transportation problem <BTP> can be stated 
as follows : a set of supplies (a ) and a set of demands <b > are 

I 1 

specified such that the total supply is equal to the total demand. 
There is a transportation time (t > associated with each supply 
point and each demand point. It is required to find a feasible 
distribution which minimizes the maximum transportation time 
associated between a supply point and a demand point such that 
distribution between two points (x^^> is positive, i.e.. 


Minimize z 


Max 1 mum C 
C ( i . j ) 


X > 0 > 

1 J 


sub ) ect 


to 

m 

2 


I I 
1 


= a 


X 


1 


1 * . • 
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m 



m n 

2 s ' 2 ^ ■ 

1=1 j = 1 


> 0 . 


Many articles on bottleneck transportation problem has 
been published in literature. A thorough review of literature has 
been provided by Szwarc C243. Hammer [2^]. Szwarc C24] and 
Srinivasan and Thompson C263 have used primal simplex approaches 
for solving BTP. They have specialized the primal simplex approac- 
hes of linear programming for the transportation problem. The 
capacitated version of the BTP has been studied by Russel et 
al . C273 . 


2.8 BOTTLENECK LINEAR PROGRAMMING PROBLEM 

Bottleneck linear programming problem <BLP> has been 
addressed by Russel et al. C273 and Frieze C283. The BLP can be 
stated as 


Minimize z = ^ x^ 

1 € C 

Subject to 

A X = b. 


c 

J 


* 

= w 


> 
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X > 0, 

where w = Maximum C c^, j : x^ > 0 }, A is an (mxn) matrix, x is 
an n vector, b is an m vector and c is an n vector. The BLP has 
wide applications in military operations, transportation of 
perishable goods C27] and assembly line balancing C283. Most 
algorithms of BLP are direct extension of the BTP. 

2.9 BALANCED OPTIMIZATION PROBLEMS 

Given a finite set E , a vector c of element costs, i.e, 

c = C c ;e e E }. and a family F of feasible subsets of E, the 
e 

balanced optimization problem is to find s which minimizes 

(maximum fc : e e S } - minimum Cc : eeS }) overall seF. 
e e 

BALANCED ASSIGNMENT PROBLEM 

The balanced assignment problem is to find a feasible 
assignment which minimizes the difference of the values of most 
expensive and least expensive costs used in the assignment. The 
balanced assignment problem has been studied by Martel lo et al. 
[293- They have specialized the feasibility questions of the 
classical assignment problem and solved balanced assignment problem 
in OCnS . 
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BALANCED LINEAR PROGRAMMING PROBLEM 

The balanced linear prograrrening problem has been studieo 
by Ahuja C503 and can be stated mathematically as follows : 

Minimize I MaximumCc x > - MinimumCc x.} 

L ] € j 1 ^ ^ ^ 

subject to 

Ax = b , 

X > 0. 

where A, b, c, x are |mxnj, |mxlj, jlxnj, jnxlj matrices respectiv- 
ely and JSN= C 1,2,. .,n >. 

Ahuja C303 has solved BLPP parametrically by transforming 
the problem into a linear program by introducing additional 
constraints and treating these constraints as parametric upper ana 
lower bounds. The algorithm is based on a relationship between BLPP 
and parametric MLPP and parametrically solves the MLPP . 



CHAPTER III 


ALGORITHMS FOR BALANCED NETWORK FLOW PROBLEM 

In this chapter we describe algorithms for the balanced 
network flow problem <BNFP>. We present two algorithms to solve 
the BNFP, namely parcanjSiTtc algorithm and biTuxry search algorithm. 
We discuss the basic ideas and concepts behind the algorithms 
before their formal statements. Later, we solve a numerical 
example to illustrate them. 


3.1 RAY SOLUTIONS 


Although one can always use linear programming approach 

to solve the BNFP, there are certain cases where due to nature- of 

the constraints, an optimum solution can be found without recourse 

to any rigorous technique. Kaplan ZZl has shown that such a case 

arises when there exists a feasible ray solxition, i.e., c^^x^^ = t 

V <i,j) € A for some scalar t i 0, Ray solution, if it exists, 

clearly solve the BNFP. To identify such a possibility, consider 

the constraints Ci . and Cl. 3^ of the BNFP . Substituting x^^ - 

t/c in these constraints, it can be easily shown that an optimum 
1 } 

ray solution of BNFP exists if and only if there exists a scalar t 
which satisfies the following two conditions; 


b< 1 ) 

I - 2 n/c,j> 

J € N j € N 


b(2) 


I - 1 

i € N ] 


(1/C -) 
) 2 

€ N 



b(n) 


y <i/c ) - y (1/C ) 

L n j L j n 

3 € N j € N 


Max 1 mum( c 

I J 

( 1 , j > <£ A 


) < t < Minimum(c u ) 
n Ml] 

( l , J ) € A 


C3. 2: 


Obviously conditions C3./J> and C3. 2J> are quite uniiKel* 
to be satisfied by most BNFP. So rigorous techniques are requirei 
to solve them. 

It IS easy to see that if the BNFP does not possess a ra 
solution, then any optimum solution is a boundary point of th 
constraints Ci.2D and C1.3S>. For, if (x*,y*, 2 *) is an optimu 
solution, which is an interior point, then there exists in th 
neighborhood of x* a solution x° defined as : 


* 

X 


^/c . 

* * 
I f C X = z 

1 J 


1 J 

13 1 J 

* 

X 

+ 

£^/C 

^ * 
i f c . X = y 

J 


1 3 

13 13 

* 

X 


a ^ 

otherwise. 

1 ) 


1 3 



Where a are such that x° is feasible to BNFP. 

1 J 

For some arbitrarily small £ > 0 , we have z - y = z 
y - 2^: < z - y , which contradicts the optimality of x . 


3.2 PARAMETRIC ALGORITHM 


Parametric algorithm of the BNFP is a spec i al i zat i on 


that for the balanced linear programming problem (BLPP> by Ahu 



C50]. The special structure of the network flow problem offers 
several benefits, particularly streamlining the computations, 
eliminating the need to explicitly maintain the simplex tableau. 
The tree structure of the basis permits the parametric algorithm to 
achieve these efficiencies. 


RELATIONSHIP BETWEEN MNFP AND BNFP 

It can be easily seen that BNFP. as already stated in 
Section 1.2. is equivalent to the following linear programming 
problem. 


Minimize <z - y) 
subject to 


C3.3:> 


wnere d 

1 1 


1/c 

I J 


'S X -Yx =b(i),V 
L 1 j L ] 1 

) eN J eN 


X < 

Minimum(d z. 

u , 

1 J 

1 J 

1 ] 

X > 

Max 1 mum (d y, 

i 

1 ] 

1 J 

1 ] 

V < 1 , J > 

e A. 



1 € N, C3.4^ 

) , V < 1 , ] ) € A , C5. 

> , V < 1 . J > € A, C3,e: 


The BNFP is closely related to the minimax networi 
flow problem MNFP(y), defined for a parameter y : 


Minimize z 
subject to 




1 J 

< 


jsN 

Minimumfd 


b< i ) , 

iJ^’ 


V 1 € N, 

> , V < 1 . J > € A , 
1 J 


C3. 7 
C3.8 
C3. 9 


C3. iO 



C3. 


X . > MaximumCd y, i ),V'Ci,])€A. 

Let x<.y) IS the optimum solution of MNFP<y> with z(.y> as 
tne objective function value. The following results establishes 
the relationship between BNFP and MNFP. 


PTpposi tion 3, i . T?i&r& exists a. y sxich that ctn optiTnxim. solxitton. of 
HNFPCyy solves BNFP. 

Proof: Let (x, y, z ) be an optimum solution of BNFP. 

can de easily seen that <x ,z >is an optimum solution of 
For if z<y*> < z* , then the solution x(y*) is feasible to to BNFP 
and has objective function value equal to z(y >“'y z -y , wnicn 

contradicts the optimality of (x ,y ,z ) for BNFP. g 


Now 1 1 
MNFP(y*> . 


Proposi.txorx 3.2. sCy2> is a pieceij>r.se linear convex fvnction in y. 

Proof: It follows from the well Known result that tne objective 
function value of linear programming problem (minimization objecti- 
ve) with parametric right hand sioe is a piecewise linear convex 
function. For details, see Murty C513. _ B 

Let 0 = e < e^ < . . . < a be the values of y where the 
1 z t 

slope of z<y> changes. Further, let the slope of z<y> m the 
interval ®k + l^ '^k each k = 1, . . - < t 1. Now the 
optimality criteria for BNFP can be stated by the following 


theorem. 
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Th&orem. 3. i . If h is t ?\9 smallest index for xjfhich ^ i, then 

optinvum. solxition of BNFP. Fxir-theT, if t/%©n 

xCe^^ solves BNFP. 

Proof: As per Proposition 5.1, we can restrict our choice to 

x(y>. BqS y < e^, while looking for an optimum solution of BNFP. 
In the BNFP, the objective function value of x<y> is M<y) = z<y)-y, 
which is a piecewise linear convex function since it is formed by 
addition of a piecewise convex function and a linear function. 
Hence a local optimum solution of M(y) is a global optimum solut- 
ion. The slope of M(y) in the interval is - 1 < 0 

and in the interval 3 i s Uj^ - 1 > 0 . Thus M(y) achieves a 

global optimum at y = e and x(e ) is an optimum solution of BNFP. 

K. JN 

Thus the first part of the theorem is proved. The second part of 
the theorem is obvious in view of the preceding discussion. H 

PARAMETRIC SOLUTION OF MNFPCY5 


Consider the following network flow problem. called 
as NFP(y,z), which treats y and z as two parameters. 


Minimize Ox 
subject to 


> X 

Z 1 ) 

=b(i), V i 

€ 

N, 


C3. 123 

jeN 

jeN 





VI 

X 

Hi n imum< d z , u ) , 

1 j i 1 

V 

{ 1 , J > 

€ A 1 

C3. 133 

X > 

1 J 

Max imum< d y , 1 > , 

U i 1 

V 

< I , J ) 

€ A. 

C3. 143 
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The above problem is a network flow problem where we are 
interested in finding a feasible flow satisfying constraints C3. iSy 
to C3.143. A basic solution to the above problem is defined by a 
triple <B.L,U>. where 8 is a spanning tree with n-1 arcs. The 

subsets B,L, and U partition the set of arcs A into three parts. 
The set L denotes the set of arcs at their lower bounds and U 
denotes the set of arcs at their upper bounds. We further divide L 
into two sets denoted by and L^, where denotes the set of 

arcs in L which are at their real lower bounds, i.e., = 1^^, 

denotes set of arcs in' L which are at their parametric lower 
bounds, i.e., x = d y. Similarly, we denote U and U as a set 

1 J I J 1 K 

of arcs in U which are at their real upper bounds and parametric 
upper bounds respectively. Now if we set 


^1 J ] 

V 

< 1 , 


- U^. 

X = d z 

1 J U 

V 

<1 . )) 

e U 


V 

( 1 , 

J) 

€ L , 
r 

T3 

it 

X 

V 

< 1 . J > 

€ L 


then there exists a unique x^^ satisfying C3.123. If this x^^ 
further satisfies C3. i 33 and C3. i 43 for all (i,]) e B, then we 
refer (B,L,U) as an optimxmt basis struct-uTS af NFPCy,z3. 


Proposition 3.3. Ci3 If HFPCy,-z3 has a feasihlo solxition^ then 
there al'UKiys exists an optiimun. basis strxictxire CBtl~,U3. 

Cii3 There always exists at least one spanning tree basis B which 
solves the HFPCy ts3 . 
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Proof: We start with x as a feasible solution to NFP(y,z). 

the flow in some arc (i,j>, x . = Minimum<d z,u >, then we put it 

1 J 1 j 1 3 

in U; else if x = Maximum{l ,d y), then we put it in L. The 

1 3 1 j 1 3 ^ ’ 

remaining arcs are put in 8. Now B may contain a number of cycles. 
We remove these cycles by selecting one such cycle and sending flow 
in clockwise <or ant i— c lockwi se > direction in it, until one arc 
reaches its lower or upper bound. Such an arc is deleted from B 
and added to either L or U depending upon its bound. We note that 
by augmenting flow in a cycle we are preserving the mass balance 
constraint at each node. By repeating the above procedure we get a 
cycle free solution from an initial feasible solution. Now if the 
basic arcs in the cycle free solution connect ail the nodes, it is 
clearly a spanning tree solution; otherwise we ado some restricted 
arcs (arcs of L and U) to the basic arcs such that the resulting 
basis connects all the nodes and it does not contain any cycle. 
Thus we a get a spanning tree basis and an optimum basis structure 
<B,L,LJ> from any feasible solution of NFP<y,2). | 


Now in the spanning t 

basis creates two subtrees. 

different subtrees form a cut. 

defines a cut. We represent a 

arc (i.j), as We denote N 

subtrees formed by removing the 

contains node i and N contains 

3 

be a set of arcs in the cut 


I 3 

head nodes 


1 I e 


i n N . 
3 


S imi lari 


ree basis, removing any arc from the 
Arcs whose nodes belong to two 
Thus each basic arc uniquely 
cut, formed by removing the basic 
and N as the set of nodes of the 


I 3 


basic arc (1,3), such that 

node 3. Clearly, N^+N^ = N. 

Q whose tail nodes lie in N 
1 3 1 

y, let Cr* be a set of arcs in 

7 . i 3 


N 

1 

Let 

and 

the 
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cut set Q whose tail nodes lie in N and head nodes lie in N . 

1 ] ) 1 

Now we can state the following relationship between the above 
defined terms. 


Q 


= ■|^< i , j > : 1 € N^and j e + ■^< i , j > : i ^ and j e 


Q = 

1 ] 1 J 1 J 


Now we establish a result which is crucial for the 
development of our algorithm. 


Theorern 3. 2. An optimxmi basis str-uctuTe CB,L,U^ of NFPCy ,z-S> with 
th& corr&s ponding solutiorx x and s^olutton x?ctlxL& 2 ? an optimxmi 

solution of HHFPCyy if and only zf thoro exists a basic arc Cp^q^ 
zn the basis B satisfying any one of the followzn^ two sets of 
condi tion: 


ci:> 

X = 

PQ 

cL y 

P<? 





C3. i 5a3 


Ci,j^ 

€ L fOT 

al 1 

<i>J^ 

€ 

pp 

C3. i 5b3 


Ci ,j:> 

« U for 

aL 1 


€ 

P9 

C3. i5c> 

c^i:> 

X = 

P<? 

d. s 





C3. i6a3 



€ L for 

all 

a 


pp 

C3. i6b3 


Ci,J^ 

G U for 

al 1 

ci,j:> 


PP 

<3. i6c3 


Proof I CHecesst tyl> We prove it by 


we s 


how that contains at ieast 


one 
* 


Let (x y be the optimum solution and z 


contradiction.- First of ail 
parametric upper bound arc. 
be the objective function 


value of MNFP(y), such that it does not have any parametric upper 



bound arc. Consider the following solution for a small value of e. 


0 * 0 * 
X = X ;2 = 2 


Clearly the solution (x°,2°> is feasible to MNFP<y) , with objective 

function value 2 - s, which contradicts the optimality of (x ,2 >. 

Hence, x contains at least one parametric upper bound arc. 

Now letd y<x <d zV<i,j>€B. Add a parametric 
I J 1 1 I J 

upper bound arc (g,h> to the basis. Adding <g,h> to the basis will 

form exactly one cycle. Now augment a flow Af = d ,£ in this 

gh 

cycle, in a direction opposite to that of (g.h), for a reasonably 
small value of e. Repeating the above steps for all the parametric 
upper bound arcs, we get a solution x° defined as below: 


0 

X 

11 

V < 1 , J ) 

€ L ; 

* 

= X 

- € d 

V 

C 1 , 

, j> 

€ U 

i J 

1 } 

1 1 

1 1 

1 J 




0 

* 

= X 

V < 1 , ] > 

« U ; 

* 

= X 


V 

< 1 , 

. J ) 

€ B 

1 ) 

1 } 

r 1 } 

1 j 

1 1 





where y = ) d ^ t , ; 

1 j L gh 1 j , gh 

(g,h)€Llp 


t = 0 if the basic arc (i,j) is not a member of the cycle 

I J .gh 

formed by adding nonbasic arc (g,h> to the basis; 


t =1 if <i,j> IS a member of the cycle and the 

ij.gh 

orientation of (i,j> is opposite to that of (g,h>; 


t = -1 if (i,]> is a member of the cycle and the orient 

iJ.gh 
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ation of (1,3) IS same as that of (g.h). 

It IS clear that the solution is feasible to tne 

MNFP<y> and has the objective function value z-e, which contradicts 
the optimality of the solution x*. Hence there must exist at least 
one basic arc, say <p,q), for which x = d y or x = d z. 

pq pq pq pq 

Let there be exactly one basic arc for which x -= d y. 

pq pq 

Further, let condition C 3 . i 5 by be violated for some arc (k,!) € 

Q . This means that (k,I> e U. Using the argument as before we 

pq 

obtain the solution x°. Consider the flow value of basic arc <p,q) 
in the solution x° : 

x° = d y + £ y where y = T d t , 
pq pq pq pq L gh pq.gh 

<g.h>€ U 

P 

If V i 0 , then x° is feasible to the MNFP(y> and has the object- 

pq 

ive function value z-e . When y ^ 0 , add the arc (k,l> to tne 

PQ 

oasis, this forms a cycle. Now augmenting a flow equals to 
in this cycle in the direction of (p,q> we obtain a solution x 
which IS defined as below: 


X 

= X . , 

V <1, 

1 J 

yv 

^kl 

1 ] 

* 

^ ^kl 

^pq 

X 

* 

= X 


i j 

1 ] 

1 J 

✓V 

X 

* 

= X 

= d y 

pq 

pq 

pq 


(L U U)-C(k. 1)>; 


y £, V<i,j)€B-C<p,q)} 
1 j ,pq ^ pq 


It can be easily verified that x is feasible to the MNFP(y) and 
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has the objective function value 2 . - s , which again contradicts 
the optimality of solution x . A similar reasoning can be used to 
show the necessity of the remaining cases. ^ 


CSxif ftciency^ Now we will show that if there exists a basic arc, 
say (p,q), in the basis which satisfies condition C3. i5J> or C3. i6^ , 
then the corresponding solution x and solution value z solves the 
MNFP(y > . 

Let a basic arc Cp,q) satisfies condition C3. i53 and the 
corresponding solution is x . Adding the mass balance constraints 
C3. 4J> for all the nodes in N we get 

q 


leN 


- X 


( 1 , J )€S~ 
pq 


pq 


- y X 

Z. 1 3 

( 1 . 3 )€(Q 


C3. 17 :> 


pq 


- C(p.q>>> 


Substituting x = d y in the above equation and rearranging both 

pq pq 


the both sides we obtain; 


leN ( 1 , 3 (i,j>€(QI^ - C(p.q)>) 

q ' pq pq 

Let x° be an optimal solution of the MNFP<y). Now we can 

o * 

state the following relationship between x and x : 


I . 


x° = X* + A. . V (i.j) € L, since x° is 
1 3 13 13 

the MNFP<y) 

x° = X - A. , V (i. 3 > € U, since x is an 
13 13 1 3 


feasible tc 


upper bouni 


2- 
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on the optimum solution of the MNFP<y> 


where > 0 V (i,j) s L u U. Similarly, adding the mass balance 

constraints C3.-4^ of all nodes in for the solution x° we get; 


pq 


= 2 


b< 1 > 
isN 


-Jcx 

( i 


+ A ) 

]>€ 

pq 


+ 


y <x* - A ) C3. iQ:> 

( 1 . J )€(Q - (p, q> > > 


On simplification, we get: 


x° = d y - y A 

pq pq L 1 j 

( 1 . ] )€ <Q - C<p,q)>) 

pq 


C3. 20:> 


It can be seen that any positive value of A^ ^ makes the 
solution x° infeasible to the MNFPCy) . Hence A =0V{i,i>€Lu 

1 j 

U. So the solution x has same objective function value as that of 

■M' 

X . A similar proof can be given when condition C3. 16^ is satisfi- 


ed . 



We call an optimum basis structure of the NFP(y,z> 
satisfying optimality condition of Theorem 5.2 as optimxmi basis 
strv.ctxire of HNFPCyJ>. We call the basic arc (p.q) satisfying the 
optimality condition as critical arc and the corresponding cut 
as critical c^t^. It may be noted that the critical cut satisfies 
the following property ; (a) All the arcs in the critical cut are 
restricted, i . e. , at their upper bounds or lower bounds. (b) All 

the upper bound arcs are in the same direction, all the lower bound 





arcs are also in the same direc’tion and it is opposite to the 

orientation of all upper bound arcs. 

Let (B.L.LI) denote the optimum basis structure of the 

MNFP<y) at y = y with Q as the critical cut. Let x be the 

- pq — 

corresponding solution and ^ be the objective function value. 
Please note that this basis structure can be obtained by solving a 
MNFP as described in Chapter 4. We now determine the maximum value 
of y, denoted by y, such that (B.L.LI) remains an optimum basis 
structure of the MNFP(y> for every y e (y, y^) . The interval (_y. y> 
IS called the characteristic interval of (B.L.U). 

DETERMINING THE CHARACTERISTIC INTERVAL 

The characteristic interval of (B.L.U) is determined from 
the following cons i derat i ons : 


(1) Increase y by an amount Ay. set x^^ = d^^(y + Ay) for each 

( 1 . j ) € L . 

P 

(2) The value of z increases by an amount, say Az . and flows on 

ail parametric upper bound arcs are set at their changing 

upper bound, i.e., x^^ = d^^(z+Az). V (i.j) € U^. 

(4) The flow value on each basic arc, say (i.j). changes uniquely 

as the y value increases by one unit. Let us denote this 

changes by s ^ ^ . 

(?) The relationship between Az and Ay is determined from the 
criteria that Q remains the critical cut. Let us denote 

pq 

the ratio Az/Ay by 6. 
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(6) The maximum permissible value of Ay is determined from the 
fact that the flow value on each basic arc remains within its 
permissible bounds. 


COMPUTING THE RATIO AZ/AY 


Consider the arcs in the critical cut. As y increases by 

an amount Ay, all the parametric lower bound arcs in the critical 

cut carry an additional amount of d Ay from one subtree to 

1 ] 

another subtree. Similarly, the parametric upper bound arcs in the 

critical cut carry an additional amount of flow d Az each in the 

1 ] 

reverse direction. In order to satisfy the mass balance constra- 
ints C3.4I> these additional incoming and outgoing flows should 
balance each other. This criteria is used to compute Az/Ay. Thus 
Az/Ay can be computed depending upon the following cases: 

Cas 9 I. WK&tx. tK& critical arc Cp,<f^ is at its paramstric Icfwsr 
bowid. 


Equating the additional incoming and outgoing flow we get 


1 


+ d Ay = y d Az 
pq Z. I J 


( , . , ) ^ n Lp 


( 1 , ] ) e n U 

pq p 


y ^ 

L I ) 


+ d 


pq 


=> At/ A y 


( 1 . J ) € Qpq n Lp 

y ^ 

< 1 , j ) S iTp n Up 


C3. 213 
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Case 1 1 . 


When the critical arc is at its parametric -upper 

bound. 


y ^ 

L 1 j 


Az/Ay 


( 1 , i ) € §“ n L 
pq p 

7 d + d 

A 1 J pq 

( 1 , j > € ^ n- u 
pq p 


C3. 223 


It may be noted that if © > 1, then as per Theorem 5.1, 
the corresponding solution x is optimum to BNFP and the algorithm 
can be terminated at this stage; otherwise it should be continued. 

Increasing y by an amount Ay causes an increase in flow 


on parametric lower 

bound arcs by d Ay 

1 J 

and 

that 

on 

parametr 1 c 

upper bound arcs by 

d Az . This increase 
I j 

1 n 

f 1 ow 

on 

parametr i c 

arcs causes excess 

or deficit at a node. 

wh 1 ch 

1 s to 

be 

balanced 


by corresponding flow changes in basic arcs. So as y increases by 
one unit, flow in each basic arc changes in a unique way. It is 
not difficult to see that as long as the current basis structure 
remains feasible, the flow on each basic (i,j> e 8 is of the form 
X + s Ay . 

-1 j 1 j 


COMPUTING FLOW CHANGES IN BASIC ARCS 


In order to determine the flow changes on a basic arc 
ti.j), as y increases by one unit, following steps can be used- 


Step 1. 


Set n. 

I 


0 V 1 e N and perform the following operations: 
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For 

each 

< I . J ) 

n 

= n - 

d 

I 

1 

I ] 

n 

ii 

d . 

j 

j 

1 J 


For 

each 

( 1 , 

. j > 

n = 

n - 

e 

d 

1 

1 


1 J 

n = 

n + 

e 

d 

J 

J 


1 J 


St&p II. In the spanning tree corresponding to the basis select a 

node of degree one, say i. and find a unique basic arc 

<g,h> such that g = i or h = i . If g = i, then set s . = 

^ ’ gh 

n ; else s , = - n . Modify n = n + n, . Delete the 

1 gh 1 g g h 

arc Cg.h) from the basis and repeat this step, until no 
arc remains in the basis. 


At y = y, flow on a basic arc, say (i,j>, reaches either 
its lower or upper bound. The arc (i,]> is then dropped from the 
basis and an eligible nonbasic arc is selected to enter the basis. 


FINDING LEAVING ARC IN THE BASIS 


A basic arc leaves the basis if it obstructs the increase 
in y value. A leaving arc Ci,)) is determined from the criteria 
that the flow on a basic arc lies within its upper and lower bound, 
i .e.. 
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Maximum<i , d (y+Ay)> < x + s Ay < MinimumCd (z+Az), u >» V 

— 1 ] i ) i j - 1 ] 

(i. }> ^ C3- 33:> 

Using the above four inequalities we obtain: 


Ay 


I j 


•1 j 


d y 
1 j ~ 


< “S 


1 J 


t d ) 
1 ] 


00 


if <-s +d > > 0 
I ] i J 

otherwi se 


Ay 


1 J 


X - d z 
1 3 - 

(-S + © d > 

1 1 I } 

00 


if (-S ta d ) < 0 
U 1 } 

otherwi se 


Ay 


1 1 




1 ~ X 

1 j -i J 


i J 


CO 


if s < 0 
1 J 

otherwise 


Ay 


I j 


u - X 
I J -i J 


I ] 


00 


if s > 0 
I J 

otherwi se 


Let Ay 


1 


I 1 


= MinimumC Ay* , Ay , Ay , Ay > 

1 J i j 1 ] 1 J 


let Ay = Ml n imumCAy ^ ^ , V (i,]> 


Further , 

B, While determining the 
permissible increase in y, i.e*. Ay, we have to consider nonbasic 
arcs as well* The Ay should be such that a nonoasic arc does not 
change from its parametric bound to its real bound or vice versa. 


for whenever such a change occurs, the s 


1 } 


values of basic arcs 


change. There are two such possibilities: 


A real lower bound arc reaches its parametric lower bound. 
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^ parametric upper bound arc reaches its real upper bound. 


Considering these possibilities we obtain: 


Ci:> 

Ay^ 

1 J 

= - y. V 

J > € L 

C2y 

Ay*^ 

1 1 

= - z)/e 

V (1. 


L 

€ U 


Let 


Ay 

Ay 

Ay 


1 


u 


Minimum(Ay! >,V(i,j)€L-L i 
1 ) P 

Minimum<Ay^^ ) , V (i.j) « ; 

Min imum(Ay^ ,Ay*^,Ay*^) . 


Now y = ^ + Ay. Clearly, <B,L,U) is an optimum basis 
structure of MNFP<y) for all y e <y.y^> with optimum objective 
function value z + 0 (y - y>. The ^ = z + 0 (y^- y> is the object- 
ive function value and x is the optimum solution at y = y. Let 
<k,l> be the arc for which minimum Ay is achieved. There are three 
possible cases: 


Ay = Ay ^ = Ay^^j^ < = > A real lower bound arc <k,i) reaches its 
parametric lower bound. 

C2J> Ay = Ay*^ = Ay!^ , < = > A parametric upper bound arc (k,l> reaches 

its real upper bound. 

C3J> Ay = Ay = Ay , < = > A. basic arc <k,l> reaches its lower or 

i\ X 

upper bound. 


If case 1 or case 2 occurs, then the following steps can 

be followed. If (k,l) e Q , i.e., an nonbas i c arc in the critical 

pq 
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cut changes its bouna, then fresh computation of the 6 and s^^ is 
required; otherwise the 9 does not changes and only updating the 
^ 1 ) are requirea. We describe below a procedure to upoate 
the s^^ values efficiently. 


Procedure 


UPDATE ; 
1 


If (Ay = Ay > then 

Let (i,]> be the arc for which Ay^ = Ay^ 

I ] 

Entering arc (i.j^ to the basis B will form a cycle 
For all basic arc (g,h> in that cycle 

If orientation of (i,j) and (g.h) are same then 


Else 


®gh ~ ^gh ^ ^i j 


s = s ~ d 
gh gh i \ 


Else 


Endif 

Let (i,j) be the arc for which 

Entering arc (i,j> to the basis B will form a cycle 
For all basic arc (g,h> in that cycle 

If orientation of (i.]> and (g.h) are same then 


®gh ~ ®gh '^i 


Else 


s,=s,+d 9 

gh gh i j 


Endif 


Endif 


If case 3 occurs, then basic arc (k.i) leaves the basis 
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depending upon the following conditions 


1 . 

“^^kl’ <k.i) leaves 

tne 

basis at 

I ts 


parametric lower bound. 





2. 

2 

If Ayj^^ = Ay^j^, then arc (k.l) leaves 

the 

basis at 

its 


parametric upper bound. 





5. 

^^kl ~ ^^kl’ 

(k , 1 > leaves the 

bas 1 s 

at Its 

real 


lower bound. 





4. 

If Ayj^j^ = Ay^j, then arc 

( k , 1 > 1 eaves the 

bas 1 s 

at its 

real 


upper bound. 






If the leaving arc 

<k , 1 ) IS the 

critical arc 

(p.q> 

itself, it implies that = 

'='ki7 = ^ 

= - z 

Thus 

the 


current solution is a ray solxition and the algorithm can be 
terminated. 

FINDING ENTERING ARC TO THE BASIS 

When the basic arc (k.l) leaves the basis, a cut is 

produced. Now a nonbasic arc is selected to enter the basis such 
that the spanning tree property of the basis is maintained and the 
entering arc does not restrict the increase in y value. An 
entering arc is selected depending upon the following cases. 

CcLs& I. Wh&n basic arc Ch,LS> l&acv&s the basis at its lower boxind. 


When the basic arc <k,l) leaves the basis at its lower 
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bound, a nonbasic arc becomes eligible to enter the basis if it is 
a member of cut any other arc enters the basis, the 

spanning tree property of the basis would not be maintained. Two 
types of arcs in the cut are eligible to enter the basis. they 

are, (1) a lower bound arc whose orientation is opposite to that of 
the arc <k,l), <2) a upper bound arc whose orientation is same as 

'that of the arc <k,l). It may be noted that if any other arc 
enters the basis it will restrict further increase in y value. So 
we define a set of eligible entering arcs K as below: 

-|^(i.j) : (i,])€Ln or ( i , ] > € U n 

Case II. When basic arc leaves the basis at its vpper bovnd. 

Similarly, we define the eligible entering arc set K as 

below: 

K= -|^(i,j) : <i,j>€Ln or < 1 , j > € U n 

While selecting an entering arc from the set K, prefer- 
ence IS given to the arcs which are not members of the critical 
cut. This ensures that the critical cut continues to satisfy the 
optimality condition after a basis change. However, when all the 
arcs in K are member of the critical cut, i.e., K ^ ^pq' then any 
arc from K is selected; in this case. though the critical arc 
remain the same, the critical cut changes. When K is empty, Q. , 
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satisfies the optimality condition, and it is made the new critical 
cut and corresponding arc is made new critical arc. 

ALGORITHM STATEMENT 


We summarize the preceding discussion by the following 
stepwise description of the algorithm. 

Step f. Set y = Mi n imumC c ^ ^ 1 ^ ^ > V (i, 3 > € A. Solve MNFP with y = 
y. If the problem is infeasible, go to step 6; else let 
X and z be the solution and the objective function value 
respectively. Obtain an optimum basis structure (B.L,U) 
with arc (p.q> as the critical arc and as the 

critical cut . 


Step 2. Compute the ratio 0 = Az/Ay depending upon the following 


cases 


Case i. Critical arc is at its parametric lower boxind. 




Pu 

< i , ) ) 


+ d 


pq 


e = 


0^ n L 
pq p 


y d 

Z 1 j . 

< 1 , j ) € n 



4 ? 


CCLS& 2. 


Critical arc is at its param&tric -app&r boxirtd, 

i • e. , X = d z 
P9 PR 


e 




< 1 , ] ) € n L 
pq p 


y 

L 1 ] 

< 1 . j > 


Q n U 
pq p 


d 

pq 


St&p 3. 


CompxLte S CFlow choTL^e in basic arc Ci , jO as y incr&as- 
&s by one for each basic arc in the following 

manner: 


CcO Set = 0 V 1 € N and compute the following: 

n = n - d and n = n + d ;V(i,j)<sL 

i i 1 j j j i ) ’ p 

n = n - 0 d and FI = n +0d ;V<i,j>€U 

i i i J J J i ] ’ P 

Cb3 In the spanning tree corresponding to the basis 

select a node of degree one. say i. and find the 

unique basic arc (g,h> such that g = i or h = i . If 

g = 1 , then set s ._ = fl ; else s , = -FI . Modify 11 

gh 1 gh i g 

= n + FI. . Delete basic arc <g,h) from the basis 

g n 

and repeat this step, until no arc remains in the 
bas 1 s . 


Step 4 Finding the leaving arc 


For each ( i , j ) e B, compute the foi lowing: 



4 a 


Ay 


1 j 




d y 
1 J 


(-S + d > 

n I ] 

00 


if (-S + d ) > 0 

1 ] 1 J 

otherwise 


Ay 


1 J 


X - d z 
~n I J 

<-s + © d ) 

1 J 1 J 

00 


If (-S. + © d > < 0 

I J 1 ] 

otherwi se 


Ay 


I 3 


= S 


i - X 
ij -1 J 


i J 


00 


1 f s <0 
1 J 

otherwi se 


Ay 


I J 


/r u - X 

I 3 -13 


1 3 


00 


if s > 0 
1 3 

otherwise 


Let Ay 
Ay 


1 3 
b 


1 Z 

Minimum (Ay ,Ay 

13 13 

= Min imum (Ay ) . 


A 5 A ^ 


( 1 . 3 ) 


B 


1 3 


For each ( 1 , 3 ) e L - L . compute Ay^ = c 1 

P 13 1313 

Let Ay^ = M 1 n imum(Ay | ^ , V (i,j> e L - L^ . 


y . 


( c u - z > 

For each ( 1 , 3 ) e U compute Ay*^ = — 

P 'J © 

Let Ay*^ = Minimum(Ay'^ > V (i,j> e U . 

13 P 


Let Ay = Min imum(Ay*^ , Ay3 , Ay^> . Set Ay , z = z + 9 

Ay, X = X + s Ay V (i,j) € B. Let (k,l> De the arc 

^ — 1 J — 1 J 13 

for which minimum Ay is achieved. 
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If Ay = Ay^ or Ay =-Ay^. then 

If Ay = Ay^, lower bound arc (k.i> reaches its 
parametric lower bouna. 

If Ay = Ay'^, then parametric upper bound arc (k,l> 

reaches its real upper bound. 

If <k , i ) € Q__. then 

pq 

go to step 2 

* 

Else use UPDATE to modify s values and 

I J 

go to the beginning of step 4. 

End i f 

Else basic arc (k,i> reaches its bound depenaing on the 
following conditions. 

If “ Ay^pthen (k,l> leaves the basis at 

parametric lower bound. 

2 

If Ay^^j = Ayj^j, then (k,i) leaves the basis at its 
parametric upper bound. 

If Ay^^^ = Ay^^ , then (k,l> leaves the basis at its 

real lower bound. 

• / A 

If Ay^^^ = '^^kl' (k,i> leaves the basis at its 

real upper bound. 

End if 

Step 5. If the leaving arc <k.l) is the critical arc <p,q> 
Itself, the current solution is a ray solution, go to 
step 7; else define an eligible entering arc set K as 


foil ows : 
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K 


" (i,j>€Ln $7 I (i.]) eUnO^, 

k 1 K 1 

: <i,j)€Ln or < i , j > € U n 

*>'ki = '^''ki “■' 


If K IS empty, make the new critical cut, go to step 
2; else find an entering arc from K. While selecting an 
entering arc, preference is given to the arcs whicn are 
not members of the critical cut. If all elements of K 
are members of the critical cut, i.e., K ^ Q . then 

pq 

select any arc from K to enter the basis. Find the new 
critical cut and go to step 2, else go to step 5 for 
fresh computations of 

Step 6. The BNFP is infeasible, STOP. 


Step 7. COptimal soLutiorO The optimal solution is = d 

( 1 . ] > € L 
L - 

B. The objective function value is z - y 


X 

P -1 J 


d^^z V <1,J> € Up, = 1^ J V (I, J> 


L,x =u V<i,j)€U-U and x V <i.j) 

p — 1 J I J P — 1 J — I j 


The preceding steps of the parametric algorithm is represented by 
the flow chart shown in Figure 5.1. 

It may be recalled, from Section 1.5, our assumptions 
that BNFP IS defined over all the arcs in a network, i.e.. A = A. 
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COMPUTE eiz/^y <#) 


arc ^ 
a MeMber of 

v crt. cutjX^ 


Critical cut 

changes, find 

new crt, cut 
__ 


^ >= i 


OptiMal sol. 
found, x*-x, 


Find 

— 4 -{ 

entr. arc Ho 


eligible^V, 
entering arc set 
ewpty ? 


COMPUTE S, , for al 1 

I V 

basic arc <i,J) 

Find the permissible 
increase i n y <^y) , 
y = y+!!iy,z = z+^iiy 


STOP 


Ray sol. obtained 

y = z 

y— — 

I 

leaving arc 
^ and crt. arc 
^\.saMe ? 


Arc set changes 


Update Sj^ 
values for 
all has. arc 


Find leaving 
Arc to basis 


Figure 3.1 Flow Chart of Parawetric ftlgorithw for BHFP 












>0 


This may not be the case in many applications as one would like to 
balance weighted arc flows in a subset A of arcs with ASA. The 
following modification while finding tne leaving arc will take care 
of such cases. Let us denote A = A - A . Set Ay^ = oo and Ay^ = 
CO for (i,j) € A . This modification will work, because the 
parametric upper bound and the parametric lower bound conditions 
will not be checked for arcs in A , whereas the real lower bound 
and upper bound conditions will restrict flows of arcs to remain 
within their respective bounds. 

3.3 BINARY SEARCH ALGORITHM 

If we represent y = Minimum( c x )V(i,j>eA and z = 

1 J I J 

Max imum< c ^ ^ X ^ ^ ) V (i.j) e A, then for a fixed value of y, say y^ , 
we can determine the corresponding minimum z value, say z^, by 
solving a MNFPCy^). Let us denote z(y) as the objective function 
value of MhJFP<y>. As per Proposition 5.2, z(y) is a piecewise 
linear convex function. The typical form of the z(y> can be 
depicted by Figure 5.2. In the BNFP, we are interested in finding 

^ ^ 'M* 

a solution x , with y = y and corresponding z(y> = z such that, 
the difference (z - y ) is minimized. The binary search algorithm 
finds y* value by using tne binary search technique in an interval 
of y. The interval is initialized as the lower and upper bound 
values of y, and is determined from the capacity constraints as 
below: 

1 .<X. <U.V<1,]>€A 

1 J 1 ) U 






Showing Rchnonihip Botw^on Objeonv* Fnnetioni'r; »f MNFP(y} andy 


Figure 3.2 


'’F 
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Figur» 3.3 FIom Chart For Binary Saaroh Algorithu 








>5 


<=> C 1 <C X <C U V(i,]>€A 
1 ) I J 1 J 1 J I J 1 J 

<=> MinimumCc 1 ) < Minimum(c x ) < Minimum<c u > 

1 J 1 J 1 J 1 ] 1 3 1 i 

<=> MinimumCc i ) < y < Minimum(c u > 

I ] I J I J I 3 

The binary search algorithm continuously improves the lower bouna 
or the upper bound in each iteration, reducing the search interval. 
Theoretically, the algorithm will never terminate, since the 
interval will not merge to a point. In order to prevent this, a 
very small interval s is to be set at the beginning. so that 
whenever the search interval becomes smaller than , the algorithm 
stops. Thus the binary search algorithm does not give an exact 
solution but gives a solution, which is within a small range ot 
tne optimal solution. 


ALGORITHM STTATEMENT 

A stepwise description of tne binary search algorithm is 
given below. 

Slept Set y^ = M i n imum( c ^ ^ 1 ^ ^ ) V (i.j> eA,y^- 
Minimum(c_^u^^> V (1.3) € A and let ^ be tne oesireo 
interval for termination. Solve MNFP(y^). If the 

problem is infeasible, then the corresponding BNFP is 
infeasible and Stop; else let x^ be the corresponding 
solution and be the objective function value. 

St.p 2 If <y^ - then go to step 5; else compute the 



^4 


f ol lowing 


= —5 + Vi 


2<y - y , ) 

u 1 


+ y 


d 5 • "1 

Solve MNFP<y„). If infeasible, then set y^ = y^, and go 


to the beginning of step 2. Solve MNFP(y^). If infeasi- 
ble, then set y^ = y^ and go to the beginning of step 2. 

Let z and z_, be the objective function values of 
c d 

tively. If “ y^^ ^ 


MNFP(y^) and MNFP(y^> respectively 

y ) then set y = y^ else set y 
' c u a 

beginning of step 2. 


Go to the 


Simp 3 y* = y, . z* = and x* = Xj. The objective function 

* 


va 


lue of BNFP IS 2 * - y* and the solution is x . Stop. 


The flow chart of the binary search algorithm is snown 1 


n 


F 1 gure 5 . 3 


3 . 4 - 


NUMERICAL EXAMPLE 


we solve a numerical example to illustrate various steps 
of the parametric algorithm. Consider the sample network shown in 
Figure 5.4. The first number besides each arc is its cost (weight) 

,-tboaritv The label associated witn 
and the second number is its capac i y. 

I / .j xsw ^ node The lower bound on 

eacn node is the supply/demand of that nooe . 

each arc is zero. 

/ t ^ - n we set y = 0. After solving 
Since Minimumlc^ ^1 j j> - 0. we sec y 
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( 6 , 6 ) 




(4,6) 


by)®- 


/5 7 ', 


;^3to) 


I l; , 


io(3X:® 


(5,7) 


'i I 


(8,6) 




(4.6) 


=^i-15 


Numerical Example Problem 


Figure 3.4 




Inirial Basis for the Example Problem 


Figure 3.5 



MNFP<y) and then transforming the solution to 
solution, we obtain a basis as shown in Figure 
structure and the objective function value *z* of 
given below: 


^6 

a spanning tree 
5 . ^ . The bas i s 
me so i ut 1 on i s 


B = C<1 . 2> , (1 .4> , (2. 3>}, L = 0, U = C < 1 , 5 ) , < 2 . 4 > . ( 5 . 4 > > , 

Up = C(1.5).<2,4)> 

Since 1^^ = 0, L = and ^ critical arc is (1,4), 

and it IS at the parametric upper bound. The critical cut Q. = 

1 4 

C<1 ,4) , (2,4) , (5,4)>, and z = 27.69. The ratio B = L'z/Ly = 0. We 

tabulate the s values of all basic arcs and the Ay values in 

M 1 ] 

Tab 1 e 5.1. 


Tctbl& 3, i 


Arc 

s 

1 ] 


Ay^ 

1 J 


Ay 

1 1 

Ay° 

(1,2) 

0 . 00 

6.00 

00 

00 

6 . 00 

6 .00 

(2,3) 

0.00 

17.84 

00 

CE> 

17.84 


(1,4) 

0 . 00 


00 

00 

27.69 



u u DU' 
Ay^,^ = Ay^^ = oo, since 0 = 0. So Ay = Minimum(Ay ,Ay ) = 6.00, y = 

y + Ay = 6.00, The basic arc (i,2) leaves the basis 'at its para- 
metric lower bound. The eligible entering arcs sec are K = 
C(l,5>>. Arc (1,3) enters to the basis. The new basis is as snown 
in Figure 3.6. The new basis structure is shown as below: 
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B = a2,5), (1,4>.(1,5)}, L = C<1.2)>, U 


C ( 2 , 4 > } . 


The recomputed s 

i ] 


and Ay 

I J 


are 


= £<2.4> . (5,4)}, 

shown in Table 5.2. 


U 

P 


Tahl& 3.2 


Arc 

s 

1 ] 

Ay' 

1 J 

. 2 

Ay 

1 ] 

Ay" 

1 ] 

Ay 

1 ] 

Ay*^ 

(2,5) 

0.00 

55.55 

14.76 

9.25 

9.25 

9.25 

(1.4) 

0.00 

21.69 

00 

CO 

21.69 


(1,5) 

0.00 

11 .85 

00 

00 

11.83 



t3 

Ay = Ay = 9.25, y = y + Ay = 15.25. The basic arc (2,5) leaves 
the basis at its upper bound. The eligible entering arc set K = <p. 
The critical cut changes now. The cut = £ ( 1 , 2 ) , ( 2 , 4 ) . ( 5 , 4 ) } . 
Any parametric arc from Q 25 ’ (1,2), enters the basis as a 
critical arc. The updated basis is shown in Figure 5.7. The new 
basis structure is: 

B = £(1 ,4) , (1 , 5) . (1 .2)}, L = .?). U = £(2,4), (5,4), (2,5)}, U = 
£(2,4)}. New 6 = 0.835. The recomputed s^^ and j j shown in 
Table 5.3. 


Table 3.3 


Arc 

s . 

1 ] 

Ay'. 

1 J 

Ay^ 

1 J 

ISH 

Ay 

i J 

Ay*^ 

(1.4) 

0.16 

5.54 

00 

00 

5.54 


(1,5) 

0 . 00 

4.77 

00 

00 

4.77 

4.77 

(1.2) 

-0.16 

t ® 

1 

74.77 

20.77 

20.77 


























\ 




Basis After 1st Pivot 

Figure 3.6 



Basis After 2nd Pivot 

Figure 3.7 
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U H 

^^24 = Minimum<Ay .Ay^> = 4.77, y = y +Ay = 20.0. z = z 

T S Ay 51.66. The basic arc (1,5) leaves the oasis at its 
parametric lower bound. The eligible entering arcs are K = 
C(2,5)}. Arc (2,5) enters to the basis, the new basis is as shown 
in Figure 5.8. The new basis structure is: 

B = Cd ,2) , (1 ,4) , (2,5>>, L = €(1.5)}. U = € ( 2 . 4) , ( 5 , 4 ) } . U 

P 

€(2,4)}. The recomputed value of & = 1 .855, so the optimal 
solution is found. The solution is tabulated in Table 5.4. 

Xcth to 3^ ^4 

z* = 51.66, y* = 20.00, z* - y* = 11.66 


Arc 

F 1 ow ( X ) 

Cost(c ) 

1 ] 

U 

1 1 

x" 

u 

(1.2) 

5.55 

6.00 

6.00 

20 . 00 

(1,5) 

4.00 

5.00 

7.00 

20.00 

(1,4) 

2 .67 

8 . 00 

5.00 

21 . 55 

(2,5) 

6 . 00 

4.00 

6.00 

24.00 

(2.4) 

6.55 

5.00 

7.00 

51 .67 

(5.4) 

1 6.00 

1 

4.00 

6.00 

24.00 































Basis After 3rd Pivot 

Figure 3.8 



CHAPTER IV 


IMPLEMENTATION DETAILS AND COMPUTATIONAL RESULTS 

In this chapter, we describe various details required in 
the implementations of the parametric algorithm, the binary search 
algorithm and the linear programming approach for balanced network 

flow problem. We also present computational results of these 
algor 1 thms . 

-4-.1 SOLVING MINIMAX NETWORK FLOW PROBLEM <MNFP) 

The parametric algorithm and the binary search algorithm 
as discussed in the previous chapter need a solution methodology to 
solve the minimax network flow problem, MNFP(y>. where y is a 
parameter. Keeping the performance and simplicity of coding in 
mind, we extended the Primal— Dual algorithm, for n^inim<xx traTtsport — 
atiorx problem, of Ahuja C25] for general networks. 

PRIMAL-DUAL ALGORITHM FOR MNFPCY3 

The MNFP<y) defined over a network G = (N.A) is to solve 
the following optimization problem. 

Minimize z 
subject to 

2 *1) ■ 2 ’‘ii ' ‘ 

j € N J € N 



62 


Maximum(l , d y> < x < Mininium<u , d z>, V 
‘ J * J 1 3 1 3 1 3 

( 1 , j ) € A , 

where = 1/c^^ for all (i,j> e A. 

First of all, the above problem is transformed into a 
maximum flow problem with parametric upper bounds by adding a 
source node s and joining it with each supply node i by an add-* 
itional arc (s,i) with capacity b(i>. Similarly, a sink node t is 
a.dded and each demand node j is connected with t by an additional 
arc <),t> with capacity -b(j>. We call these additional arcs as 

t i i 

artificioil. arcs . Let G = (N ,A ) be the enlarged network. Now 
consider the following parametric maximum flow problem with z as 
a parameter. 


Maximize v 


sub 3 ect to 

3 € N 3 € N 


V 1 f 1 = s 

- 0 1 f 1 ^ s , t 

-V if 1 = t 


Maximumdj^, ^ 


< Minimum<u , d z> 
1 } 13 


( 1 . 3 > e A , 

0 < X ^ b<i) V 1 € N and b<i) > 0, 
s 1 

0 < x.^ ^ -b<j) V 3 € N apd. b( 3 > < 0, 


V 


The above problem is a maximum flow problem with para 
metric capacity. It is well known that the objective function 
value of a linear programming problem (maximization ob 3 ective> with 
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a 1 gor 1 thm 


Step i 


Step 2 


Step 3. 


Step 4 


IS given below. 


Compute g 


b( 1 ) 

y cl . 

L 1 j 
j € N 


for 1 € N such that b< i ) > 0 


K - -b<]> 

h^ - fQP J ^ [vj gych that b< ] ) < 0 


2^. 


) € N 


Set k = 0, z = Maximum Maximum(g^), Maximum(h^) ^ 


1 : leN and ] : jeN and 

b( 1 ) >0 b( ] > <0 


k k 

Define capacity of each arc <i,j) s A as cl^^z if ci^^z < 

u else u . Similarly, define lower bound of each arc 
1 J 1 J 

as the larger among (d^^y, Maximum flow 

# 

algorithm to determine the maximum flow in G from source 

k k 

to sink. Let x be the resulting flow and v be the 

1 J 


max imum 


flow value. If v = T. go to step 4. 


Label those nodes in G for which flow can be augmented 

. * 

from node s. Let C = C <i,]) = (i,j) € A , i is labeled 

k 

and j IS not labeled >. Let S = C (i,j> : (itj) -e C n A 

and d < u }. Compute w*^ = T d and + (T 

1 J 1 ) Z- 1 J 

( 1 , J >€S 

- v'^>/w*^. k = k + 1 and go to step 1. 

# k 

(Optimal solution found) z = z, x = x and Stop. 

The above steps of the algorithm is represented by flow 
4.2. The primal-dual algorithm needs a 


chart as shown in Figure 



START 







COHPUTE 

T = Zb<i) ; b(i) > 0 
gj= b(i)/Edj^,i : b<:i)>0 
hj= b<J)/Zdjj j ; b(J><0 

k = 0 

z* = HaX<MdX, ( 9 ,) 

Define lower bound of <I,J) 
AS MaxiwuMCdj jy, 1 , j)for all 
arcs 


♦ 



Flo« Cb.rt ot Pri«al-Dual 


Algorithm iot MHFPCsJ 


Figure 4-2 
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maximum flow subroutine to solve the MNFP<y), and its computational 
performance is determined by the efficiency of the maximum flow 
code. We use a version of maximum flow subroutine wnicn is cased 
on shortest axigmenting path algorithm. 


4-.2 GETTING THE INITIAL BASIS AND CRITICAL CUT 


The primal-dual algorithm, when terminates, does not 
yield a basic solution (extreme point solution). This is in fact 
true for most network optimization algorithms because most network 
flow problems do not possess a unique optimum solution. Rather, 
they have alternate optimum solutions. We need to convert the 
optimum solution obtained from primal-dual algorithm to an extreme 
point solution, i.e., a cycle free solution in the network optimiz- 
ation language. We use the following steps to obtain a cycle free 
solution. 


Step i After solving the MNFP(y), add all arcs at their upper 
bounds (either parametric or real) to set U. Add all 
arcs at their lower bounds (parametric or real) to set L. 
and the remaining arcs to set B. 


Step 2 


If B does not contain a 
identify a cycle in set B 
B by augmenting flow in 
direction until one arc in 
bound (parametric or real) 


cycle, go to step 5; else 
and then remove the cycle from 
clockwise (or ant i -c 1 ocKw i se ) 
the cycle reaches its upper 
or its lower bound (parametric 
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or real). Delete the arc from B and add it to U or L 
depending on its bound. Repeat this step. 


Step 3 A cycle free solution has been found. Stop. 


At the end of the above steps we get a cycle free 
solution, i.e., the set B consists of arcs which form a forest 
(union of trees). The trees in B are disconnected by cuts. We now 
require to identify a critical cut among these cuts and aod some 
restricted arcs between these cuts so as to get a basis whicn is a 
spanning tree. We use the following steps to get a spanning tree 
basis and identify a critical cut. The following procedure 
partitions the set N of node into two parts: marked nocas and 
unmarked nodes. 


Step 1, 


Make all nodes in N unmarked. Select any node from N 
say 1 , and mark i t . 


Step 2 Use search aleorithn to identify ail nodes which are 
reachable from i. Make ail these nodes marked. 


Step 3 


If there is no 
cr 1 1 1 cal cut i s 
ity condition 
unmarked) . I f 


unmarked node, then go to step 5. 

yet to be found, then check the optimai- 
of Theorem 5.2 for the cut C(marked, 
It satisfies, then it is the critical cut. 


Step 4 


If CCmarked. unmarked) is the critical cut. 


then add a 
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parametric arc ; else add any arc from this cut to the 
basis. In the former case the arc is called critical 
arc. Go to step 2. 

Step 5 A critical cut and a spanning tree basis has been 
obtained. Stop. 

We use breadth first search version of the search 
algorithm to perform various operations like identifying a cycle, 
finding a cut. The basic iterative steps of the breadth first 
search procedure is described as follows. 

Algorithm breadth first search; 
be^irx. 

unmark all nodes in N; 
select any node i ; 
mark node i; pred( i ) = 0; 

Queue = C i > ; 

Whtle Queue do 

begtn 

select first element of the queue; 

if node 1 IS incident to arc (i,j> and ] pred( i ) then 

begin 

mark node j ; 
pred ( ) > = 1 ; 
add node j to Queue; 


0nd; 
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olso dGlGte nod© i "from Qusug; 


©nd; 


&nd; 


4-3 


NETWORK REPRESENTATION 


Network representation is one important factor that 
significantly affects the running time of an algorithm. A clever 
network representation not only reduces the storage requirement, 
but also facilitates easy updating, and easy maintenance of the 
intermediate results. We use the forward star reoresentat i on to 
store networks as it suits the requirements of our algorithm. We 

describe below this representation. 

In the forward star representation, we arrange the arcs 
in A in certain order; we first arrange arcs emanating from node 1. 
then from node 2, and so on. Arcs emanating from tne same node can 
Be numoered arbitrarily. We then sequent i a 1 1 y stc.-e their tail, 
head, cost and capacity in that order. We maintain a pointer with 
each node l. denoted by pointd). that indicates the smallest 
number of the arc emanating from node i. Hence the outgoing arcs 
from node i are stored at position pointci) to point<itl)-l m the 

arc list. Ifpoint<i)> point(.fl> - K t*'- ‘ 

mnintCi) = i and point(n4*l) 
t 1 s t ency » we set pointvi/ i 

outgoing arc. For cons i sT.ei>>-y . 

, 4 . ,ro thp Table 4.1 forward star representat i or 

= m+1. We Illustrate in the laoie 

of network already shown in Figure 5.4. 
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Table 4. i 


Arc no. 

Point 

Arc 

Cost 

Capac 1 ty 

1 

1 

(1,2) 

6.0 

6.0 

2 

4 

(1,3) 

3.0 

7.0 

5 

6 

<1 ,4) 

4.0 

6.0 

4 

7 

(2.5) 

4.0 

6.0 

5 

7 

(2,4) 

3.0 

7.0 

6 


(3.4) 

4.0 

6.0 


4-A 


DATA STRUCTURE USED 


The parametric algorithm maintains a basis at every 
Iteration. While performing various operations on the basis, like 
Identifying a cycle, finding a cut. we require access to all arcs 

incident to a node. We use the following data structure to 

maintain the basis and perform above operations efficiently. We 
use an array, called bstar. which arranges the arc numoers of arcs 
in the basis in certain order. It first arranges all arcs incident 

to node 1. then node 2. and so on. It stores same arcs in two 

places because it is incident to two nodes. Similar to the forward 

star representation, we keep a pointer to each node in a array. 

^+..,..■1- ,r>n nnsition of arcs emanating from 
called bpoint . which gives starting position o 

^ 4- 1 n-p the basis, we keep ali basic arcs in 
a node. For easy updating of 

w and keep track of the leaving arc from a 

another array, called bas . and k p 

called bdeiet, whose i-th 
basis by maintaining a logical array. 

^ , leaves the basis. We modify the 

element is made true if arc i la 
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bstar and bpoint array whenever the basis changes. 

^•5 IMPLEMENTATION OF LINEAR PROGRAMMING CODE 

We compared the execution times of parametric and binary 
search algorithms with the simplex algorithm applied to the 
enlarged linear programming problem. We used the XMP library 
developed by Roy E. Marsten C52] for solving the enlarged linear 
program. The XMP library contains subroutines written in VANILLA- 
FORTRAN to solve different sizes of linear programs. We used those 
subroutines which are based on the prirna.1 simplex algoT-ttlxm ana 
exploit the sparseness of the constraint matrix. The number of 
eligible nonbas i c variables which are put in the candidate list and 
the re factor 1 zat 1 on frequency of the basis which have an important 
bearing on the performance of the code were kept 6 and 50 respecti- 
vely. The above values were set so that the LP code gave a better 
performance . 

4.6 CHOICE OF NETWORK GENERATOR 

For large scale computational testing, we need a network 
generator to generate random networks which are possioly good 
representation of reality. We used Netgen, the network generator 
program developed by Klingman, Napier and Stutz in 1974 C533. 

Netgen can generate different sizes of minimum cost flow, transpor- 

tation and assignment networks. 

We generated the arc capac.t.ee, costs, total supply .n 
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the following manner. The capacity of an arc was generated as a 
random number in the interval C0.5n, 2. On], cost in the interval 
Cl.O, m] , and the total supply is set equal to 50n, where n anc m 
are number of nodes and number of arcs in the network respectively. 
We set the total number of source node equal to 0.2n, the total 
number of sink node equal to 0.25n. The total number of transship- 
ment source nodes and transshipment sink nodes were set half of the 
total source and total sink nodes respectively. The percentage of 
skeleton arcs which were given maximum cost was 40% and the 
percentage of capacitated arcs was 40%. 


A-.7 COMPUTATIONAL RESULTS 


The parametric algorithm, and the binary search aigontnm 
were coded in FORTRAN-77. We also used the FORTRAN codes for the 
primal simplex algorithm of XMP library. The three cooes were 
extensively tested on HP 9000/850s computer systems running unoer 
HP-UX operating system. We tested the codes with different sizes 
of problems generated by netgen. We cons.aereo seven sizes of 
problem and for eacn problem size we solved five prooiems. We 

noted down the CPU time and the following counts for each code. 

Far pcrrcmj»trtc algorithm : 

1> Time to get initial basis. 

2) Time for parameterization, 

5) Number of pivot operations. 


For binary saarch algorithm : 
1> Number of search points 
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2) Nurftbsx^ of C3.iis to ths rn3.x i rnurn “flow rout ins 
5) Total number of flow augmentations 
4) Total time to obtain the solution 
For prinval simplex algorithm : 

1) Total number of pivot operations 

2) Total time to obtain the solution 

We have tabulated these values, averaged over five 
problems and presented in Tables 4.2 to 4.4. All time information 
are given in seconds. For comparison purposes we have given the 
CPU times of all the three codes in Table 4.5. For better insight 
of the data, the CPU time verses problem size plot for all tne 
three codes have been shown in Figure 4.5 to 4.5. The average 
number of iterations verses problem size plot for parametric 
algorithm is shown in Figure 4.6. Using the regression analysis, 
we have estimated the time and number of iterations the parametric 
code will take for different problem sizes. The ratio of estimated 
time and actual time, and estimated iterations and actual number of 
Iterations for different problem sizes have oeen shown in Figures 


4.7 to 4.12. 
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TablB 4. 2 


Computational Results for Parametric Algorithm 


No . of 

Nodes 

No . of 

Arcs 

No . of 

Herat ions 

Minimum 

T ime 

Maximum 

T ime 

r 

Time to 

Get basis 

Time for Para- 
meter i zat 1 on 

Total 

T ime 

10 

20 

1.20 

0.01 

0.05 

0.016 

0.040 

0.02 

10 

50 

5.20 

0.04 

0.06 

0,028 

0.022 

0.05 


72 

5.00 

0.05 

0.12 

0.052 

0.028 

0.06 

25 

50 

1 .4 

0.05 

0.07 

0,048 

0.012 

0.06 

m 

125 

5.4 

0.10 

0.27 

0.072 

0.078 

0.15 

251 



250 

4,6 

0.15 

0.51 

0.166 

0.154 

0.50 

50 

100 

5.0 

0.15 

0.51 

0.152 

0.053 

0.19 

50 

250 

12.8 

0.47 

1.05 

0.280 

0.460 

0.74 

50 

500 

12.2 

0.95 

2.04 

0.558 

0.792 

1 . 55 

100 

200 

5.0 

0.44 

0.72 

0.408 

0.102 

0.51 

100 

500 

6.2 

0.69 

1.71 

0.710 

0.450 

1 .14 

100 

1000 

22.8 

5.86 

5.56 

1.664 

2.846 

4.51 

200 

400 

5.6 

1 . 58 

2.17 

1.418 

0,262 

i . 68 

200 

1000 

51.4 

5.05 

12.19 

5.408 

4.552 

7.76 

200 

2000! 

1 

22.8 

6.25 

25.64 

5.998 ! 

5.962 

11 .96 

500 

1000 

6.4 

7.56 

10.26 

7.508 

1 . 182 

8.49 

B 

2000 

72.0 

20.58 

44.49 

15.506 

22.514 

57.82 

500 


110.8 

40.56 

140.07 

42.182 

71.658 

115-84 

1000 

2000 

79.6 

26.05 

78.59 

25.476 

50.924 

>6.40 

1000 

4000 

146.6 

67.61 

189.45 

65.256 

92.674 

155.91 

1000 

10000 

248.6 

467.54 

529.16 

164.482 

550.618 

495.01 























































































































Table 4.3 


Ih 


Computational Results for Binary Search Algorithm 


No . of 

No , of 

Ml n imum 

Maximum 

Average 

No*of sea- 

No. of call 

No. of flow 

Nodes 

Arcs 

T ime 

T ime 

T ime 

rch points 

to max. flow 

Augmentat i ons 

10 

20 

0.06 

0.65 

0.26 

B 

95.4 

735.4 

10 

50 

0.58 

0.87 

0.71 

42.4 

150.8 

1914.8 

10 

72 

1.01 

1 .41 

1.12 

166.6 

3465.0 

5465.0 

25 

50 

0.63 

1.85 

1.22 

52.0 

172.4 

3891 .2 

25 

125 

2.04 

2.91 

2.50 

46.2 

212.6 

7646.4 

25 

250 

5.10 

5.02 

5.83 

45.6 

210.0 

0172,6 

50 

100 

1 .65 

5.07 

2.99 

41 .8 

225.8 

9871 .8 

50 

250 

6.08 

7.68 

6.72 

49,8 

297.2 

8775.2 

50 

500 

8.52 

12.89 

10.49 

46.6 

261.4 

8256.4 

100 

200 


13.19 

6.75 

29.8 

189.4 

18124.2 

100 



20.97 

t 6.57 

52.2 

337.6 

6139.3 

1 

100 

1000 


34.82 

27.60 

52.8 

562.0 

4046.8 

200 

400 

9.87 

16.42 

12.99 

51.2 



219.2 

57216.4 


1000 

25.77 


36.18 

52.6 

395.0 

100610.4 

200 

2000 

52.82 


60.10 

51.2 

599.8 

148511 .0 







1 


500 

1000 

32.43 

58.04 

44.27 

55.8 

283.4 

114264.0 

^00 

2000 

74.18 

119.17 

99.34 

48.60 

425.8 

58175 . 6 

500 



218.12 

171.24 

47.4 

426.4 

577655.3 

1000 

2000 

150.88 

164.42 

146.00 

47.2 

421.8 








576.4 

511092 . ** 

1000 

4000 

188.58 

220.57 

204.49 

51.0 







504.2 

845066.2 

1000 

10000 

355.50 

466 . 1 1 

401.70 

51.6 


















































































































Table 4. 4 


Computational Results of Primal Simplex Algorithm 


No. of No. of Minimum Maximum Average No. of 

Nodes Arcs Time Time Time Pivots 



20 


50 


72 

50 

125 

250 

100 

250 

500 


200 


500 

1000 

400 

1000 

2000 


0.27 


1.28 


2.57 
1.83 
6.37 
12.21 
6.97 
20.94 
88 .53 


10.90 


30 .33 


117.78 

25.98 

74.03 

442.73 


0.40 


1.80 


3.17 

2.49 

10.03 

27.39 

7.88 

30.80 

101.37 


0.34 


1 .49 


2.89 

2.04 

7.83 

20.61 

7.33 

26.90 

94.95 


31.80 


81.20 



14 

12 

.91 

252. 

50 

68 

31 

25 

250 . 

00* 

44 

153 

.11 

601 . 

00 

. 32 

27 

.28 

250, 

. 00* 

.58 

76 

.22 

250 

. 00*! 

.93 

524.46 

1015 

.00 


Mote ♦ Opt,e,al eolut.pn could not found due to 
large number of degenerate pivots. 



























Tahl& 4.5 


Comparison of Exocuiion Times 



No .of I No . of 
Nodes 


10 


10 


10 


25 


25 


25 


50 


50 


50 


100 


100 


100 


200 


200 


200 


500 

500 

500 

1000 

1000 

1000 


20 


50 


72 


50 


25 


250 


100 


250 


500 


200 


500 


1000 


400 


1000 


2000 


1000 

2000 

5000 

2000 

4000 

10000 


Average Time 


0.26 


0.71 


1 .12 


1 .22 


2.50 


5.85 


2.99 


6.72 


10.49 


6.75 


16.57 


27.60 


12.99 


56.18 


60.10 


44.27 
99.54 
171 . 24 
146.00 
204.49 
401.70 


Parametric 

0. 

02 

0 . 

05 

0 . 

06 

0 . 

06 

0 . 

15 

0. 

50 

0 . 

19 

0 . 

74 

1 . 

55 

0 . 

51 

1 . 

14 

4. 

51 

1 . 

68 

7. 

.76 


.96 


8.49 

57.82 

115.84 

56.40 

155.91 

495.01 


0.54 


1.49 


2.89 


2.04 


7.85 


20.61 


7.55 


26.90 


94.95 


12.91 


5 


155.11 


27.28* 


76.22* 


524.46 


Note * Optimal solution could not found due to large 
number of degenerate pivots. 
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0.1424*n**0.738 
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The data clearly shows the superiority of the parametric 

algorithm over binary search and primal simplex algorithms. The 

average number of iterations for the parametric algorithm is 

considerably smaller than that of the primal simplex algorithm. 

The parametric algorithm also requires less storage requirement 

than the primal simplex algorithm. While parametric algorithm 

requires k^m + k^n storage, the simplex algorithm requires k^m^ + 
2 

k^mn + k^n storage, where k^ . ..,k^ are some constants. We have 
also noted during testing that the primal simplex code could not 
get optimal solution for many problems due to the large number of 
degenerate steps. In fact, when the problem size became larger 
than 200 nodes and 2000 arcs, the simplex algorithm could not find 
the optimal solutions for most problems. For smaller size problem- 
s, the binary search and simplex code took large time compared to 
the parametric code. For large size problems the binary search 
code became competitive with the parametric code for higher m/n 
ratios. This is because of the following reason: as the problem 
size increases, the number of pivot operations also increases for 
parametric algorithm whereas the number of search points for the 
binary search algorithm does not change so much since it depends 
upon the maximum cost and maximum capacity of the arc set. we note 
that the binary search algorithm does not give an exact solution. 
Further, the parametric code can be made faster than our implement- 
ation by choosing better data structure, whereas the binary search 
code is more or less optimized. 
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4.8 CONCLUSIONS 

In this thesis we have described two algorithms for tns 
balanced network flow problem. We have specialized the parametric 
algorithm of balanced linear programming problem of AhujaC30] for 
network problems. The parametric algorithm exploits the special 
structure of the constraint matrix to simplify the computations, 
easy updating and maintenance of the basis resulting in lesser CPU 
time. The binary search algorithm uses binary search technique to 
get near to the optimal solution. It does not give an exact 
solution but gives a solution which is within a specifiea interval 
around the optimal solution. The primal simplex algorithm to solve 
the balanced network flow problem is found to be uneconomical from 
execution time as well as storage point of view. Our computational 
testing shows the superiority of parametric algorithm over the 


other two methods. 
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